详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

内装式空射运载火箭重力出舱机箭耦合动力学分析

研究了以运输机为平台的内装式运载火箭空射过程载机和火箭的耦合动力学建模。建模针对两个阶段：第一阶段，火箭固定于载机机舱内，两者构成一个整体，按照普通刚体的力学方法处理；第二阶段，火箭解锁后，沿着舱内的发射筒向外滑行，载机和火箭形成两刚体相互作用的耦合系统，基于牛顿-欧拉法建立系统动力学模型。载机在空射火箭过程中，油门和升降舵满偏，在加速前飞的同时拉大姿态俯仰角，火箭在自身重力分量和惯性力的作用下，沿着机舱内的发射筒加速向外滑行，直至与载机分离。数值仿真分析了空射过程载机的重要力学参数的变化过程，验证了载机操控策略的可行性和安全性，可为未来中国空射运载火箭技术设计提供数据参考。     

Effect of bristle pack position on the rotordynamic characteristics of brush-labyrinth seals at various operating conditions

Over the last few decades, the research on the effect of bristle pack position on the rotordynamic characteristics of the brush-labyrinth seals is not sufficient. To this end, two kinds of brush-labyrinth seals for the bristle pack element installed upstream of the labyrinth teeth named BSU and installed downstream of the labyrinth teeth called BSD were used to investigate the effect of bristle pack position on the rotordynamic characteristics of the brush-labyrinth seals. Using the numerical model combining the porous medium model and the whirling rotor method, the rotordynamic characteristics of the BSU and BSD at various operating conditions including four kinds of pressure ratios, five kinds of inlet preswirl speeds and four kinds of rotor spinning speeds were conducted. The obtained results show that the effects of operating conditions on rotordynamic coefficients for the different seal configurations are different. The direct stiffness, cross-coupled stiffness and direct damping of the BSU are lower than those of the BSD. The rotordynamic coefficients of the BSU are more insensitive to the operating conditions variation. From the perspective of the seal stability, the BSU is a better brush-labyrinth seal configuration at high pressure ratio, high positive preswirl or high rotor spinning speed conditions. While in the case of low pressure ratio, low positive preswirl or low rotor spinning speed conditions, the BSD is a better choice.     

基于CFD高精度算法的雷电电磁场研究

雷电空间电磁场一直是电磁领域研究的热门，成熟的算法例如有限时域差分法（FDTD），传输线矩阵法（TLM）以及频域的矩量法（MoM）在计算雷电问题方面都有广泛的应用。由于计算流体力学（CFD）中的Euler方程与电磁学中的Maxwell方程有着相同的守恒形式，而且采用间断伽辽金方法（DG）已经在流场问题上得到广泛的尝试，因此引入了基于计算流体力学的DG方法来离散时域Maxwell方程，并采用网格分区并行技术加速计算，使用基于DG的圆球雷达散射截面积（RCS）算例进行测试，数值结果一致表明DG算法在求解电磁场问题上的可行性，之后通过计算一段近场雷电通道的电场分布并与解析解、某算法仿真解对比，数据基本吻合，说明该方法适合于雷电电磁场的计算。     

A progressive approach to predict shot peening process parameters for forming integral panel of Al7050-T7451

In this paper, a progressive approach to predict the multiple shot peening process parameters for complex integral panel is proposed. Firstly, the invariable parameters in the forming process including shot size, mass flow, peening distance and peening angle are determined according to the empirical and machine type. Then, the optimal value of air pressure for the whole shot peening is selected by the experimental data. Finally, the feeding speed for every shot peening path is predicted by regression equation. The integral panel part with thickness from 2 mm to 5 mm and curvature radius from 3200 mm to 16000 mm is taken as a research object, and four experiments are conducted. In order to design specimens for acquiring the forming data, one experiment is conducted to compare the curvature radius of the plate and stringer-structural specimens, which were peened along the middle of the two stringers. The most striking finding of this experiment is that the outer shape error range is below 3.9%, so the plate specimens can be used in predicting feeding speed of the integral panel. The second experiment is performed and results show that when the coverage reaches the limit of 80%, the minimum feeding speed is 50 mm/s. By this feeding speed, the forming curvature radius of the specimens with different thickness from the third experiment is measured and compared with the research object, and the optimal air pressure is 0.15 MPa. Then, the plate specimens with thickness from 2 mm to 5 mm are peened in the fourth experiment, and the measured curvature radius data are used to calculate the feeding speed of different shot peening path by regressive analysis method. The algorithm is validated by forming a test part and the average deviation is 0.496 mm. It is shown that the approach can realize the forming of the integral panel precisely.     

星载激光通信终端在轨机动对温度影响

对激光通信终端在轨瞬态温度变化开展了仿真，以期研究在轨机动的影响、热分析和热试验时机动模式的简化模拟。通过合理地分析与简化，建立终端的轨道热分析模型，准确模拟在轨机动，根据典型机动模式、外热流和涂层退化等因素设计了计算工况，得出了终端在轨运行过程中的温度场随时间和姿态变化规律。结果表明：不同机动模式下的温度变化存在差异，最大可达23.0℃，采用固定姿态或一维转动的简化热分析或热试验不能准确模拟实际飞行温度，甚至不能部分替代二维转动热分析或热试验；光学天线和反射镜的温度控制是制约终端工作的瓶颈，为终端设计合适的避光机动策略，可大幅度提高温度稳定度和均匀度。研究结果可以为光机电设备在轨机动策略的设计和改进提供参考。     

机载托架自动校准系统的设计与实现

为了解决机载托架传统校准方法精度低、耗时、费力等问题,设计了一种高精度、高效率并且操作简单的成品托架自动校准系统。该系统运用调平精度较高的"循环多次"最高点不动调平方法,建立高精度机载托架校准的静力学数学模型,运用VC与Matlab的COM接口编程技术开发机载托架自动校准系统,并建立良好的人机交互界面。试验结果表明,利用该方法对机载托架进行校准,无论是校准精度还是校准效率都较传统方法有很大提高。     

Fault mode probability factor based fault-tolerant control for dissimilar redundant actuation system

This paper presents a Fault Mode Probability Factor (FMPF) based Fault-Tolerant Control (FTC) strategy for multiple faults of Dissimilar Redundant Actuation System (DRAS) composed of Hydraulic Actuator (HA) and Electro-Hydrostatic Actuator (EHA). The long-term service and severe working conditions can result in multiple gradual faults which can ultimately degrade the system performance, resulting in the system model drift into the fault state characterized with parameter uncertainty. The paper proposes to address this problem by using the historical statistics of the multiple gradual faults and the proposed FMPF to amend the system model with parameter uncertainty. To balance the system model precision and computation time, a Moving Window (MW) method is used to determine the applied historical statistics. The FMPF based FTC strategy is developed for the amended system model where the system estimation and Linear Quadratic Regulator (LQR) are updated at the end of system sampling period. The simulations of DRAS system subjected to multiple faults have been performed and the results indicate the effectiveness of the proposed approach.